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ABSTRACT 

We examine the relationship between the frequency of close companions (stellar and planetary 
companions with orbital periods < 5 years) and the metallicity of their Sun-like (~ FGK) hosts. We 
confirm and quantify a ~ 4er positive correlation between host metallicity and planetary companions. 
We find little or no dependence on spectral type or distance in this correlation. In contrast to the 
metallicity dependence of planetary companions, stellar companions tend to be more abundant around 
low metallicity hosts. At the ~ 2a level we find an anti-correlation between host metallicity and the 
presence of a stellar companion. Upon dividing our sample into FG and K sub-samples, we find a 
negligible anti-correlation in the FG sub-sample and a ~ 3a anti-correlation in the K sub-sample. A 
kinematic analysis suggests that this anti-correlation is produced by a combination of low-metallicity, 
high-binarity thick disk stars and higher-metallicity, lower-binarity thin disk stars. 

Subject headings: binaries: close - stars: abundances - stars: kinematics 



1. INTRODUCTION 

With the detection to date of more than 160 exoplanets 
using the Doppler technique, the observation of Gonza- 
lez (1997) that giant close-orbiting exoplanets have host 
stars with relatively high stellar metallicity compared to 
the average field star has gotten stronger (Reid 2002; 
Santos et al. 2004; Fischer & Valenti 2005; Bond et al. 
2006). To understand the nature of this correlation be- 
tween high host metallicity and the presence of Doppler- 
detectable exoplanets, we investigate whether this corre- 
lation extends to stellar mass companions. 

There has been a widely held view that metal-poor stel- 
lar populations possess few stellar companions (Batten 
1973; Latham et al. 1988; Latham 2004). This may have 
been largely due to the difficulty of finding binary stars in 
the galactic halo, e.g. Gunn & Griffin (1979). Duquen- 
noy & Mayor (1991) investigated the properties of stellar 
companions amongst Sun-like stars but did not report a 
relationship between stellar companions and host metal- 
licity. Latham et al. (2002) and Carney et al. (2005) 
reported a lower binarity for stars on retrograde Galac- 
tic orbits compared to stars on prograde Galactic orbits 
but found no dependence between binarity and metallic- 
ity within those two kinematic groups. Dall et al. (2005) 
speculated that the frequency of host stars with stellar 
companions may be correlated with metallicity in the 
same way that host stars with planets are. 

In this paper we describe and characterize the correla- 
tion between host metallicity and the fraction of plane- 
tary and stellar companions. In Section 2 we define our 
sample of close planetary and stellar companions and we 
describe the variety of techniques used to obtain metal- 
licities of stars that do not have spectroscopic metallic- 
ities from Doppler searches. In Section 3 we analyze 
the distribution of planetary and stellar companions as 
a function of host metallicity. We confirm and quantify 
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the correlation between planet-hosts and high metallicity 
and we find a new anti-correlation between the frequency 
of stellar companions and high metallicity. In Section 4 
we compare our stellar companion results to analogous 
analyses of the Nordstrom et al. (2004) and Carney et al. 
(2005) samples. 

2. THE SAMPLE 

We analyze the distribution of the metallicities of FGK 
main-sequence stars with close companions (period < 5 
years). For this we use the sample of stars analyzed by 
Grether & Lineweaver (2006). This subset of 'Sun-like' 
stars in the Hipparcos catalog, is defined by 0.5 < B — 
V < 1.0 and 5A(B -V) + 2.0 < M v < 5.4(5 -V)- 0.5. 
This forms a parallelogram -0.5 and 2.0 mag, below and 
above an average main-sequence in the HR diagram. The 
stars range in spectral type from approximately F7 to 
K3 and in absolute magnitude in V band from 2.2 to 7.4. 
From this we define a more complete closer (d < 25 pc) 
sample of stars and an independent more distant (25 < 
d < 50 pc) sample. See Grether & Lineweaver (2006) for 
additional details about the sample definition. 

2.1. Measuring Stellar Metallicity 

The metallicity of most of the extrasolar planet hosts 
have been determined spectroscopically. We analyse the 
metallicity data from three of these groups: (1) the Mc- 
Donald observatory (hereafter, McD) group (e.g. Gon- 
zalez et al. 2001; Laws et al. 2003), (2) the European 
Southern Observatory (hereafter, ESO) group (e.g. San- 
tos et al. 2004, 2005), and (3) the Keck, Lick and Anglo- 
Australian observatory (hereafter, KLA) group (Fischer 
& Valenti 2005; Valenti & Fischer 2005). 

All three of these groups find similar metallicities for 
the extrasolar planet target stars that they have all ob- 
served as shown by the comparisons in Fig. 1. Apart 
from the ~ 1000 KLA target stars analyzed with a consis- 
tently high precision by Valenti & Fischer (2005), many 
nearby (d < 50 pc) FGK stars lack precise metallicities if 
they have any published measurement at all. A smaller 
sample of precise spectroscopic metallicities has also been 
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Fig. 1. — Exoplanet Target Stars Metallicity Comparison. We 
compare the spectroscopic exoplanet target metallicities of the 
McD, ESO and KLA groups. The 59 red dots compare the ESO to 
the McD values of exoplanet target metallicity that these groups 
have in common. We find that the ESO values are on average 
0.01 dex smaller than the McD values with a dispersion of 0.05 
dex. Similarly the 99 green dots compare the KLA values to the 
average 0.01 dex smaller ESO values with a dispersion of 0.06 dex. 
The 56 blue dots compare the KLA to the average 0.01 dex larger 
McD values with a dispersion of 0.06 dex. A solid black line shows 
the slope-one line with dashed lines at ±0.1 dex. The three lin- 
ear best-fits for these three comparisons are nearly identical to the 
slope-one line and almost all scatter is contained within 0.1 dex. 
The relationship between the McD and ESO values is very close 
with a marginally looser relationship to the KLA values. Thus, 
these values for exoplanet target metallicity are consistent at the 
~ 0.1 dex level. 

published by the ESO group for non-planet hosting stars 
(Santos et al. 2005). 

Since the large sample of KLA stars has been taken 
from exoplanet target lists it also has the same biases. 
This includes selection effects (i) against high stellar 
chromospheric activity (ii) towards more metal-rich stars 
that have a greater probability of being a planetary host 
and (iii) against most stars with known close [9 < 2") 
stellar companions. We need to correct for or minimize 
these biases to determine quantitatively not only how the 
planetary distribution varies with host metallicity but 
also how the close stellar companion distribution varies 
with host metallicity, that is, we need metallicities of 
all stars in our sample in order to compare companion- 
hosting stars to non-companion-hosting stars, and to 
compare the metallicities of planet -hosting stars to the 
metallicities of stellar-companion-hosting stars. 

In addition to the metallicities reported by the McD, 
ESO and KLA groups, we use a variety of other sources 
and techniques to determine stellar metallicity although 
with somewhat less precision. These include other 
sources of spectroscopic metallicities such as the Cayrel 
de Strobel et al. (2001) (hereafter, CdS) catalog, metal- 
licities derived from uvby narrow-band photometry or 
broad-band photometry and metallicities derived from 
a star's position in the HR diagram. The precision of 
the spectroscopic metallicity values in the CdS catalog 
are not well quantified. However, many of the stars in 
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Fig. 2. — Metallicity values from exoplanet spectroscopy com- 
pared to four other methods of obtaining stellar metallicities. We 
compare the exoplanet target spectroscopic metallicities (plotted 
on the x axis as a 'reference') with the following test samples plot- 
ted on the y axis: (1) CdS spectroscopic metallicities (red dots), (2) 
uvby photometric metallicities (green dots) (3) broad-band pho- 
tometric metallicities (blue dots) and (4) HR diagram K dwarf 
metallicities (aqua dots). The mean differences between the test 
and the reference sample metallicities ([Fe/H]tcst — [Fc/H] ro f) are 
—0.05, —0.08, 0.01, and —0.10 dex respectively, with dispersions of 
0.08, 0.11, 0.14 and 0.14 dex respectively. Comparing these mean 
differences and dispersions we find that the mean differences are 
within la of the solid black slope-one line and thus we regard the 
systematic offsets as marginal. The four linear best-fits for these 
four comparisons (shown by the four colored lines) do not show 
significant deviation from the slope-one line (black) except for the 
metallicities derived using broad- band photometry (dark blue line). 

the catalog have several independent metallicity values 
which we average, excluding obvious outliers. To de- 
rive metallicities from uvby narrow-band photometry we 
apply the calibration of Martell & Smith (2004) to the 
Hauck & Mermilliod (1998) catalog. We also use val- 
ues of metallicity derived from broad-band photometry 
in (Ammons et al. 2006). For stars with 5.5 < M v < 7.3 
(K dwarfs) the relationship between stellar luminosity 
and metallicity is very tight (Kotoneva et al. 2002). Us- 
ing this relationship, we derive metallicities for some K 
dwarfs from their position in the HR diagram. 

To quantify the precision of their metallicities, we com- 
pare in Fig. 2 the different methods of determining 
metallicity. We use the high precision exoplanet target 
spectroscopic metallicities from the McD, ESO and KLA 
surveys (or the average when a star has two or more val- 
ues) as the reference sample. We compare these metal- 
licities with metallicities of the following test samples: 
(1) CdS spectroscopic metallicities, (2) uvby photomet- 
ric metallicities, (3) broad-band photometric metallici- 
ties and (4) HR Diagram K dwarf metallicities. 

The result of this comparison is that the uncertainties 
associated with the high quality exoplanet target spec- 
troscopic metallicities of McD, ESO and KLA groups 
are the smallest, with the CdS spectroscopic metallicities 
only slightly more uncertain. The uncertainties associ- 
ated with the uvby photometric metallicities are inter- 
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TABLE 1 

Stellar Samples Used in Our Analysis 









Range 






Stars with [Fe/H] Measurements 




Sample 


B 


- V 


(PC) 


Total a 


Binary 


' 1 'arfpts c 


Planet Hosts' 3 


[Fe/H] Source 


Our FGK 


0.5 


- 1.0 


d < 25 


453 


45 (9.9%) 


379 (84%) 


19 (5.0%) 


Mostly Spec. 




0.5 


- 1.0 


25 < d < 50 


2745 


107 (3.9%) 


1597 (58%) 


36 (2.3%) 


Mostly Phot. f 


Our FG 


0.5- 


- 0.75 


d < 25 


257 


27 (10.5%) 


228 (89%) 


13 (5.7%) 


Mostly Spec. g 




0.5 - 


- 0.75 


25 < d < 50 


1762 


76 (4.3%) 


1167 (66%) 


32 (2.7%) 


Mostly Phot. h 


Our K 


0.75 


- 1.0 


d < 25 


196 


18 (9.2%) 


151 (77%) 


6 (4.0%) 


Mostly Spec. 1 




0.75 


- 1.0 


25 < d < 50 


983 


31 (3.2%) 


430 (44%) 


4 (0.9%) 


Mostly Phot J 


GC k FGK 


0.3- 


- 1.0 


d < 40 


1375 


378 (27.5%) 






uvby Phot. 




0.5 - 


- 1.0 


d < 40 


1289 


346 (26.8%) 






uvby Phot. 


GC FG 


0.5- 


- 0.75 


d < 40 


1117 


291 (26.1%) 






uvby Phot. 


GC K 


0.75 


- 1.0 


d < 40 


172 


55 (32.0%) 






uvby Phot. 


CL 1 AFGK 


0.0 


- 1.0 




963 


254 (26.4%) 






Spec. 



Note. — For notes e - j, "HP Spec." is high precision exoplanet target spectroscopy, "CdS Spec." is Cayrel de Strobel et al. (2001) 
spectroscopy, u uvby Phot." is uvby photometry, "BB Phot." is broad-band photometry and "HR K Dwarf is the method for obtaining 
mctallicities for K dwarfs from their position in the HR diagram (Kotoncva et al. 2002). 

a Total number of Hipparcos sun-like stars ("H" in Fig. 4). 

b Subset of total stars that are hosts to stellar companions ("S" in Fig. 4). The percentages given correspond to the fraction S/H. 
c Subsct of total stars that are exoplanet target stars ("T" in Fig. 4). The percentages given correspond to the fraction T/H. 
d Subsct of target stars that are exoplanet hosts ("P" in Fig. 4). The percentages given correspond to the fraction P/T. 
c 63% HP Spec, 12% CdS Spec, 20% uvby Phot., 1% BB Phot, and 4% HR K Dwarf. 
f 19% HP Spec, 5% CdS Spec, 55% uvby Phot., 17% BB Phot, and 4% HR K Dwarf. 
S63% HP Spec, 18% CdS Spec, 19% uvby Phot. 

h 26% HP Spec, 7% CdS Spec, 61% uvby Phot., 6% BB Phot, and < 1% HR K Dwarf. 
! 64% HP Spec, 5% CdS Spec, 20% uvby Phot., 3% BB Phot, and 8% HR K Dwarf. 
>6% HP Spec, 3% CdS Spec, 44% uvby Phot., 36% BB Phot, and 11% HR K Dwarf. 

k "GC" is the Geneva-Copenhagen survey of the Solar neighbourhood sample (Nordstrom et al. 2004). We only include those binaries 
observed by CORAVEL between 2 and 10 times (see Section 4). 

'"CL" is the Carney-Latham survey of proper-motion stars (Carney et al. 2005). We only include those stars on prograde Galactic orbits 
(V > —220 km/s). The CL sample also includes 231 stars from the sample of Ryan (1989). 



mediate with broad-band photometric and HR diagram 
K dwarf mctallicities being the least certain. 

2.2. Selection Effects and Completeness 

To minimize the scatter in the measurement of stellar 
mctallicity while including as many stars in our samples 
as possible, we choose the metallicity source from one of 
the five groups based upon minimal dispersion. Thus we 
primarily use the spectroscopic exoplanet target metal- 
licities. If no such value for metallicity is available for 
a star in our sample we use a spectroscopic value taken 
from the CdS catalog, followed by a uvby photometric 
value, a broad-band value and lastly a HR diagram K 
dwarf value for the metallicity. We use the dispersions 
discussed above as estimates of the uncertainties of the 
metallicity measurements. 

Almost all (453/464 = 98%) of the close (d < 25 
pc) sample and 2745/2832 = 97% of the more distant 
(25 < d < 50 pc) sample thus have a value for metallic- 
ity. Given the different uncertainties associated with the 
five sources of metallicity, the close sample has more pre- 
cise values of metallicity than the more distant sample. 
The dispersion for the close and far sample are 0.07 and 
0.10 dex respectively. See Table 1 for details. 

We also investigate the color or host mass dependence 
of the host metallicity distributions. Thus we split our 
close and far samples which are defined by 0.5 < B — V < 
1.0 into 2 groups, those with 0.5 < B - V < 0.75 which 
we call FG dwarfs and those with 0.75 < B — V < 1.0 
which we call K dwarfs. This split is shown in Table 1. 
In this table we also show the total number of stars in 
the sample that have a known value of mctallicity and 
the fraction that are close binaries, exoplanet target stars 



and exoplanet hosts. 

In order to determine whether there is a real physical 
correlation between the presence of stellar or planetary 
companions and host mctallicity we need to show that 
there are only negligible selection effects associated with 
the detection and measurement of these two quantities 
that could cause a spurious correlation. In Section 2.3 
we show that the planetary companion fraction should 
be complete for planetary companions with periods less 
than 5 years for the sample of target stars that are be- 
ing monitored for exoplanets. This completeness helps 
assure minimal spurious correlation between the proba- 
bility of detecting planetary companions and host metal- 
licity. 

The stellar companion sample is made up of two sub- 
samples: those companions detected as part of an exo- 
planet survey and those that were not. The target list 
for exoplanets is biased against stellar binarity as dis- 
cussed in Grether & Lineweaver (2006). We show that 
there is negligible bias between the probability of detect- 
ing stellar companions and host metallicity in two ways: 
(1) by showing that our close sample of stellar compan- 
ions is nearly complete and (2) by using the Geneva- 
Copenhagen survey (hereafter, GC) of the Solar neigh- 
bourhood (Nordstrom et al. 2004) sample of stars, con- 
taining similar types of stars to those found in our sam- 
ple, as an independent check on our results. 

The GC sample of stars which contains F0-K3 stars 
is expected to be complete for stars with stellar com- 
panions closer than d < 40 pc. For the close sample of 
stars (d < 25 pc), the northern hemisphere of stars with 
close stellar companions is approximately complete. The 
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TABLE 2 

Metallicity and Frequency of Hosts with Close Planetary and 
Stellar Companions 



Companions 


Range 


Total 


Metal-poor 


Sun-like 


Metal-rich 


Planets 


d < 25 


19 


2 (11%) 


5 (26%) 


12 (63%) 


Stars 


d < 25 


45 a 


25 (56%) 


14 (31%) 


6 (13%) 


Planets 


25 < d < 50 


36 


3 (8%) 


9 (25%) 


24 (67%) 


Stars 


25 < d < 50 


107 b 


55 (51%) 


36 (34%) 


16 (15%) 



a An additional 2 hosts with unknown metallicity have stellar companions. 
b An additional 3 hosts with unknown metallicity have stellar companions. 



southern hemisphere of stars is also nearly complete if we 
include the binary stars from Jones et al. (2002) that are 
likely to fall within our sample (Grether & Lineweaver 
2006). We then find that ~ 10% of stars have stellar com- 
panions with periods shorter than 5 years. If we make 
a small asymmetry correction (to account for the south- 
ern hemisphere not being as well monitored for binaries) 
we find that <~ 11 ± 3% of stars have stellar companions 
within this period range (Grether & Lineweaver 2006). 
We also compare our sample with that of the "Carney- 
Latham" survey (hereafter, CL) of proper-motion stars 
(Carney & Latham 1987; Carney et al. 1994) in Section 4. 
The CL sample also contains ~ 11% of stars with stellar 
companions with periods shorter than 5 years (Latham et 
al. 2002) . We tabulate the properties of all these samples 
in Table 1. 

2.3. Close Companions 

The close companions included in our d < 25 pc and 
25 < d < 50 pc samples are enclosed in a rectangle of 
mass-period space shown in Fig. 3. These companions 
have primarily been detected using the Doppler tech- 
nique but the stellar companions have been detected with 
a variety of techniques not exclusively from high preci- 
sion exoplanet Doppler surveys. Thus we need to con- 
sider the selection effects of the Doppler method in order 
to define a less-biased sample of companions (Lineweaver 
& Grether 2003). Given a fixed number of targets, the 
"Detected" region should contain all companions that 
will be found for this region of mass-period space. The 
"Being Detected" region should contain some but not 
all companions that will be found in this region and the 
"Not Detected" region contains no companions since the 
current Doppler surveys are either not sensitive enough 
or have not been observing for a long enough duration 
to detect companions in this regime. Thus as a conse- 
quence of the exoplanet surveys' limited monitoring du- 
ration and sensitivity for our sample we only select those 
companions with an orbital period P < 5 years and mass 
M 2 > 0.001M o . 

In Grether & Lineweaver (2006) we found that com- 
panions with a minimum mass in the brown dwarf mass 
regime were likely to be low mass stellar companions seen 
face on, thus producing a very dry brown dwarf desert. 
We also included the 14 stellar companions from Jones 
et al. (2002) that have no published orbital solutions but 
are assumed to orbit within periods of 5 years. We find 
one new planet and no new stars in our less biased rect- 
angle when compared with the data used in Grether & 
Lineweaver (2006). This new planet HD 20782 (HIP 
15527) (indicated by a vertical line through the point 
in Fig. 3), has been monitored for well over 5 years but 
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Fig. 3. — Masses and periods of close companions to stellar hosts 
of FGK spectral type. We split the close companion sample into 
3 groups defined by the metallicity of their host star: metal-poor 
([Fe/H] < -0.1), Sun-like (-0.1 < [Fc/H] < 0.1) and metal-rich 
([Fc/H] > 0.1) which arc plotted as white, grey and black dots 
respectively. The larger points are companions orbiting stars in 
the more complete d < 25 pc sample, while the smaller points are 
companions to stars at distances between 25 < d < 50 pc. We di- 
vide the stellar companions into those not monitored by one of the 
exoplanet search programs (shown with an 'X' behind the point) 
and those that are monitored. Both groups of stellar companions 
are distributed over the entire less-biased region (enclosed by thick 
line). Hence any missing stellar companions should be randomly 
distributed. For multiple companion systems, we select the most 
massive companion in our less-biased sample to represent the sys- 
tem. 



only has a period of ~ 1.6 years and a minimum mass of 
1.8Mj up placing it just between the "Detected" and "Be- 
ing Detected" regions. While most planets are detected 
within a time frame comparable to the period, the time 
needed to detect this planet was much longer than its 
period because of its unusually high eccentricity of 0.92 
(Jones et al. 2006). We thus have two groups of close 
companions to analyse as a function of host metallicity - 
giant planets and stars. 

In Fig. 3 we split the close companion sample into 
3 groups defined by the metallicity of their host star: 
metal-poor ([Fe/H] < -0.1), Sun-like (-0.1 < [Fe/H] < 
0.1) and metal-rich ([Fe/H] > 0.1) which are plotted as 
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Fig. 4. — Metallicities ([Fe/H]) of 453 stars in our close d < 25 pc 
sample. Top: The lightest shade of grey are the Hipparcos 'Sun- 
like' stars in our close sample. The next darker shade of grey are 
all of the stars that are exoplanet targets. This is followed by a 
still darker shade of grey which are hosts of stellar companions. 
The darkest shade of grey are exoplanet hosts. In the Bottom plot, 
the fraction of target stars, stellar companion hosts and planetary 
companion hosts are shown by squares, triangles and circles re- 
spectively. The linear best-fit to the target fraction is shown by a 
dotted line. The linear and exponential best-fits to the stellar and 
planetary companion fractions are shown by dashed and solid lines 
respectively. 

white, grey and black dots respectively. Fig. 3 sug- 
gests that the hosts of planetary companions are gen- 
erally metal-rich whereas the hosts of stellar companions 
are generally metal-poor. Table 2 and Fig. 4 confirm the 
correlation between exoplanets and high-metallicity and 
indicate an anti-correlation between stellar companions 
and high metallicity. 

3. CLOSE COMPANION - HOST METALLICITY 
CORRELATION 

We examine the distribution of close companions as a 
function of stellar host metallicity in our two samples. 
We do this quantitatively by fitting power-law and ex- 
ponential best-fits to the metallicity data expressed both 
linearly and logarithmically. We define the logarithmic 
[Fe/H] and linear Z/Z metallicity as follows: 



[Fe/H] = log(Fe/H) - log(Fe/H) Q = log(Z/Z Q ) (1) 

where Fe and H are the number of iron and hydrogen 
atoms respectively and Z = Fc/H. We examine the 
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Fig. 5. — Same as Fig. 4 except that metallicity is plotted 
linearly as Z/Zq. All of the metal-rich (Z/Zq > 1.8) sample stars 
are being monitored for exoplanets but as the stellar metallicity 
decreases so does the fraction being monitored. This is because of 
a bias towards selecting more metal-rich target stars for observation 
due to an increased probability of planetary companions orbiting 
metal-rich host stars. The linear best-fit to the target fraction is 
shown by a dotted line. The linear and power-law best-fits to the 
stellar and planetary companion fractions are shown by dashed and 
solid lines respectively. 

close planetary companion probability P p i a net and the 
close stellar companion probability -P s tar as a function of 
[Fe/H] in Figs. 4 and 6 for the d < 25 and 25 < d < 50 pc 
samples respectively. Similarly we also examine P p i an et 
and P st ar as a function of Z/Z Q in Fig. 5 and 7, which 
is effectively just a re-binning of the data in Fig. 4 and 
Fig. 6. We then find the linear best-fits to the planetary 
and stellar companion fraction distributions as shown by 
the dashed lines in Figs. 4-7. 

We also fit an exponential to the [Fe/H] planetary (as 
in Fischer & Valenti 2005) and stellar companion fraction 
distributions in Figs. 4 and 6 and equivalently a power- 
law to the data points for the Z/Z Q plots, Figs. 5 and 7. 
The two linear parameterizations that we fit to the data 
are: 



^i„Fc/H = a[Fc/H]+P 
Pii„z/z o =A(Z/Z ) + (P -A) 
and the two non-linear parameterizations are: 

P EP = P 1O Q [ Fc / H l 



(2) 
(3) 



(4) 
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Fig. 6. — Same as Fig. 4 except for the 2745 stars in the more 
distant 25 < d < 50 pc sample. It is harder to detect distant 
planets because of signal to noise considerations which limit ob- 
servations to the brighter stars. This fainter, more distant sample 
relies more on photometric metallicity determinations than does 
the closer, brighter sample which has predominantly spectroscopic 
metallicity determinations (see Table 1). The fraction of stars be- 
ing monitored for cxoplancts is much lower than in Fig. 4. 
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Fig. 7. — Same as Fig. 6 except that metallicity is plotted 
linearly as Z/Zq analogous to Fig. 5. In this more distant sample 
we find the same trends as in Fig. 5 but they are not as prominent. 

= P (Z/Z )« (5) 

where Pq is the fraction of stars of solar metallicity (i.e. 
[Fe/H] = and Z/Z Q = 1) with companions. 

If the fits for the parameters a, A and a are consistent 
with zero then there is no correlation between the frac- 
tion of stars with companions and metallicity. On the 
other hand, a non-zero value, several sigma away from 
zero suggests a significant correlation (a, A or a > 0) or 
anti-correlation (a, A or a < 0). 

The best- fit parameters a, A and Pq (but not a) de- 
pend upon the period range and completeness of the 
sample. In order to compare the slopes from different 
samples, we parametrize this dependence in terms of the 
average companion fraction -P avg for the sample, i.e., if 
the average companion fraction for a sample is twice as 
large as for another sample, the best-fit slopes a and A 
as well as the fraction of stars of solar metallicity Pq 
will also be twice as large. To compare samples with dif- 
ferent Pavg, we scale the best-fit Eqs. 2-5 to a common 
average companion fraction by dividing each equation 



by -Pavg- Thus, we scale the best-fit parameters a, A and 
Pq by dividing each by P avg . These scaled parameters 
are then referred to as a' = a/P avg , A' = A/P avg and 
Pq = P Q /F avg . We list the unsealed best-fit parame- 
ters a, A, Pq, a along with P avg for each sample in Table 
3. The parameters a' — a/P avg and a of the different 
samples are compared in Fig. 10. 

We consistently find in Figs. 4 and 6 that metal-rich 
stars are being monitored more extensively for exoplan- 
ets than metal-poor stars as quantified by the "Target 
Stars" / "Hipparcos Sun-like Stars" ratio. This is be- 
cause of a bias towards selecting more metal-rich stars 
for observation due to an increased probability of plan- 
etary companions orbiting metal-rich host stars. Note 
that this bias is well-represented by a linear trend as 
shown by the dotted best-fit line in these figures, and is 
not just a case of a few high metallicity stars being added 
to the highest metallicity bins. We correct for this bias 
by calculating "P/T" not "P/H" for each metallicity bin. 

We find a correlation between [Fe/H] and the presence 
of planetary companions in Fig. 4. The linear best-fit 
(Eq. 2) has a gradient of a = 0.18 ± 0.07 (x? cd = 1-21) 
and thus the correlation is significant at the 2a level. The 
non-linear best-fit (Eq. 4) is a = 2.09±0.54 (x r 2 cd = 0.16) 
and thus the correlation is significant at slightly more 
than the 3a level. 

Similarly we find a correlation between linear metal- 
licity Z/Zq and the presence of planetary companions in 
the same data re- binned in Fig. 5. The linear best-fit 
(Eq. 3) has a gradient of A = 0.07 ± 0.03 (x 2 od = 1.25) 
and the non-linear best-fit (Eq. 5) has an exponent of 
a = 2.22 ± 0.39 (x 2 c d = 1-00) which are non-zero at the 
~ 2cr and ~ ha significance levels respectively. These 
results are summarized in Table 3. We can compare the 
non- linear best-fit (Eq. 5) for linear metallicity Z/Z 
and the non-linear best-fit (Eq. 4) for log metallicity 
[Fe/H] since both contain the parameter a. As shown by 
the x P er degree of freedom x?ed' the- non- linear good- 
ness of fit is better than the linear goodness of fit. We 
rely on the best fitting functional form which is the non- 
linear parameterization of our results although we use 
both parameterizations in our analysis. 

We combine these two non-independent, non-linear 
best-fit estimates by computing their weighted average. 
We assign an error to this average by adding in quadra- 
ture (1) the difference between the two estimates and (2) 
the nominal error on the average. Thus our best esti- 
mate is a — 2.2±0.5. Hence the correlation between the 
presence of planetary companions and host metallicity is 
significant at the ~ 4cr level for a non-linear best-fit and 
at the ~ 2a level with a lower goodness-of-fit for a linear 
best-fit in our close, most complete sample. 

In Fig. 4, and its re-binned equivalent Fig. 5, we 
find an anti-correlation between the presence of stel- 
lar companions and host metallicity. The linear stellar 
companion best-fits have gradients of a = —0.14 ± 0.06 
( x 2 Qd = 3 00 ) and A = _o.06 ±0.03 (x? od = 0.91) respec- 
tively, both significant at the ~ 2a level. The non-linear 
best-fit to the stellar companions as a function of [Fe/H] 
in Fig. 4 is a = -0.86 ± 0.10 (x 2 od = 1.33) and the non- 
linear best-fit to the stellar companions as a function of 
Z/Z Q in Fig. 5 is a = -0.47 ± 0.18 (x 2 od = 0.40). Av- 
eraging these two as above we obtain, —0.8 ± 0.4, which 
is significant at the ~ 2a level. All these best-fits are 
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Fig. 8. — Same as Fig. 4 for the stars in our close d < 25 pc sam- 
ple but only for FG dwarfs (B — V < 0.75). All stars have known 
metallicity in this sample. There is no apparent anti-correlation 
between metallicity and the presence of stellar companions. 
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Fig. 9. — Same as Fig. 4 for the stars in our close d < 25 pc 
sample but only for K dwarfs (B — V > 0.75). This plot shows 
a strong anti-correlation between metallicity and the presence of 
stellar companions. 



summarized in Table 3. 

Having found a correlation for planetary companions 
and an anti-correlation for stellar companions in our close 
sample and having found them to be robust to different 
metallicity binnings, we perform various other checks to 
confirm their reality. We check the robustness of both 
results to (i) distance and (ii) spectral type (~ mass) of 
the host star. 

To check if these anti-correlations have a distance de- 
pendence, we repeat this analysis for the less complete 
25 < d < 50 pc sample. As shown by the best-fits in 
Figs. 6 and 7 and summarized in Table 3 we find only 
a marginal anti-correlation between the presence of stel- 
lar companions and host metallicity for the linear best- 
fits. The non- linear best-fits however still suggest an 
anti-correlation with a = —0.59 ± 0.12 for log metal- 
licity [Fe/H] and a = -0.44 ± 0.12 for linear metallicity 
Z/Zq, which are significant at the 4a and 3a levels re- 
spectively. Combining these two estimates as described 
above we find a = — 0.5 ± 0.2 significant at the 2a level 
in the 25 < d < 50 pc sample. 

The correlation between the presence of planetary 
companions and host metallicity for the less complete 
25 < d < 50 pc is significant at the 4cr and 3a levels for 
the linear best-fits a and A respectively. The non-linear 
best-fit correlation has a = 2.56 ± 0.45 for log metallicity 
[Fe/H] and a = 3.00 ± 0.46 for linear metallicity Z/Z 



which arc significant at the 5a and 6a levels respectively. 
Combining these two estimates we find the weighted av- 
erage as above of a = 2.8±0.6 significant at the 4a 
level. 

Having found the correlation for planetary companions 
and the anti-correlation for stellar companions robust to 
binning but less robust in the less-complete more distant 
sample, we test for spectral type (~ host mass) depen- 
dence. We split our sample into bluer and redder subsam- 
ples to investigate the effect of spectral type on the close- 
companion/host-mctallicity relationship. We define the 
bluer subsample by B— V < 0.75 (FG spectral type stars) 
and the redder subsample by B — V > 0.75 (K spectral 
type stars). Since B — V has a metallicity dependence, 
a cut in B — V will not be a true mass cut, but a di- 
agonal cut in mass vs metallicity. Thus, interpreting a 
B — V cut as a pure cut in mass introduces a spurious 
anti-correlation between mass and metallicity. 

The linear best-fit to the stellar companions of the FG 
sample (d < 25 pc) has a normalised gradient of a' = 
(-0.01 ± 0.08)/10.5% = -0.1 ± 0.8 and the non-linear 
best-fit is a = -0.2 ± 0.4 as shown in Fig. 8. Both of 
these best-fits are consistent with the frequency of stellar 
companions being independent of host metallicity. The 
linear best-fit to the stellar companions of the K sample 
has a gradient of a' = (-0.27±0.07)/9.2% = 2.9±0.8 and 
the non-linear best-fit is a = —1.0 ±0.1 as shown in Fig. 
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TABLE 3 

Best-fit Trends for Close-Companion Host-Metallicity Correlation (see Fig. 10) 
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d < 25 


[Fe/H] 


Fig. 4 


Planets 


0.18 ±0.07 


4.8 ±1.2 


2.09 ±0.54 


4.5 ± 1.3 


5.0 




d < 25 


[Fe/H] 


Fig. 4 


Stars 


-0.14 ±0.06 


7.3 ± 1.4 


-0.86 ±0.10 


7.8 ± 1.3 


9.9 




d < 25 


Z/Z 


Fig. 5 


Planets 


0.07 ±0.03 


3.8 ±2.5 


2.22 ±0.39 


5.3 ± 1.4 


5.0 




d<25 


z/z 


Fig. 5 


Stars 


-0.06 ± 0.03 


8.8 ±3.5 


-0.47 ±0.18 


8.8 ± 1.5 


9.9 




d < 25 


Avg. 


- 


Planets 


- 


- 


2.2 ±0.5 


- 


- 
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d<25 


Avg. 


- 
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- 


- 


-0.8 ±0.4 


- 


- 


FGK 


25 < d < 50 


Fe/H 


Fig. 6 


Planets 


0.04 ±0.01 


1.6 ±0.4 


2.56 ± 0.45 


2.3 ±0.5 


2.3 




25 < d < 50 


Fe/H 


Fig. 6 


Stars 


-0.00 ±0.02 


3.3 ±0.4 


-0.59 ±0.12 


3.5 ±0.4 


3.9 




25 < d < 50 


Z/Zq 


Fig. 7 


Planets 


0.03 ±0.01 


2.0 ±0.6 


3.00 ± 0.46 


2.0 ±0.4 


2.3 




25 < d < 50 


Z/Z Q 


Fig. 7 


Stars 


-0.01 ± 0.01 


3.4 ±0.8 


-0.44 ±0.12 


3.6 ±0.4 


3.9 




25 < d < 50 


Avg. 


— 


Planets 


— 


— 


2.8 ±0.6 


— 


— 




25 < d < 50 


Avg. 


— 


Stars 






—0.5 ± 0.2 








d < 25 


Fe/H 


Fig. 8 


Planets 


n oo _l_ n no 


6.3 ± 1.7 


2.66 ± 0.62 


5.0 ± 1.7 


5.7 


Our 


d < 25 


TP ^ 1 U 

r ej Ti 


Fig. 8 


Stars 


-0.01 ±0.08 


10.3 ±2.0 


-0.18 ±0.37 


10.7 ±2.1 


1U.0 


FG 


25 < d < 50 


Fe/H 




Planets 


0.05 ±0.02 


2.0 ±0.5 


2.63 ± 0.46 


2.8 ±0.6 


2.7 




25 < d < 50 


Fe/H 




Stars 


0.04 ±0.02 


4.4 ±0.6 


-0.09 ±0.18 


4.7 ±0.6 


4.3 




d < 25 


Fe/H 


Fig. 9 


Planets 


0.11 ±0.10 


3.6 ±2.0 


1.57 ±0.85 


5.1 ±3.1 


4.0 


Our 


d < 25 


Fe/H 


Fig. 9 


Stars 


-0.27 ±0.07 


5.7 ± 1.9 


-0.95 ±0.14 


7.1 ± 2.3 


9.2 


K 


25 < d < 50 


Fe/H 




Planets 


0.03 ±0.03 


0.9 ±0.6 


2.10 ±2.16 


1.3 ± 1.0 


0.9 




25 < d < 50 


Fe/H 




Stars 


-0.06 ± 0.02 


2.0 ±0.5 


-1.12 ±0.19 


2.0 ±0.5 


3.2 


GC C FGK 


d < 40 


Fe/H 


Fig. 12 


Stars 


-0.15 ± 0.05 


24.0 ±1.5 


-0.28 ±0.07 


23.9 ±1.4 


26.8 


GC FG 


d < 40 


Fe/H 


Fig. 13 


Stars 


-0.05 ± 0.06 


25.1 ± 1.6 


-0.08 ±0.10 


25.5 ± 1.6 


26.1 


GC K 


d < 40 


Fe/H 


Fig. 14 


Stars 


-0.46 ± 0.10 


22.3 ±3.6 


-0.52 ±0.11 


23.7 ±3.9 


32.0 


CL d AFGK 




Fe/H 




Stars 


-0.07 ±0.06 


21.3 ±3.4 


-0.12 ±0.09 


21.6 ±3.1 


25.7 


Combined" 


Fe/H 


Fig. 15 


Stars 


-0.10 ± 0.03 


20.5 ± 1.3 


-0.22 ±0.05 


20.7 ± 1.2 


25.7 



a "Type" refers to whether the data is binned in [Fe/H] or in Z/Zq. For data binned in [Fc/H] the linear slope is a and for data binned 
in Z/Zq the linear slope is A (see Eqs. 2-5). 

b "Pavg" is defined as the number of stars with stellar companions divided by the total number of stars. We use this parameter to scale 
a, A and Pq which are then referred to by a' = a/P avg , A' = A/P avg and Pg = Pq/P &vs and can be compared between different samples 
(see text Section 3). 

C "GC" is the Geneva-Copenhagen survey of the Solar neighbourhood sample (Nordstrom et al. 2004). We only include those binaries 
observed by CORAVEL between 2 and 10 times (see Section 4). 

d "CL" is the Carney-Latham survey of proper-motion stars (Carney et al. 2005). We only include those stars on prograde Galactic orbits 
(V > —220 km/s) and with [Fe/H] > —1.3. This also includes stars from the sample of Ryan (1989). 

"Combined sample of (i) our (d < 25 pc F7-K3, Fig. 4) sample, (ii) the volume-limited GC (d < 40 pc F7-K3, Fig. 12) sample and (iii) 
the prograde Galactic orbits from Carney et al. (2005). Fig. 15 shows all three datasets. 



9 for the close d < 25 pc stars. Both of these best-fits 
show an anti-correlation between the presence of stellar 
companions and host metallicity at above the 3a level. 
Less significant results are obtained for the 25 < d < 50 
FG and K spectral type samples with stellar companions 
as shown in Table 3. 

These results suggest that the observed anti- 
correlation between close binarity and host metallicity 
is either (i) real and stronger for K spectral type stars 
than for FG stars or (ii) due to a spectral-type dependent 
selection effect. 

Under the hypothesis that the anti-correlation between 
host metallicity and binarity is real for K dwarfs, there is 
a possible selection effect limited to F and G stars that 
could explain why we do not see the anti-correlation as 
strongly in them. Doppler broadening of the line profile, 
due to random thermal motion in the stellar atmosphere 
and stellar rotation, both increase in more massive F 
and G stars due to their higher effective temperature 
and faster rotation speeds compared with less massive K 
stars. This wider line profile for F and G stars results in 
fewer observable shifting lines thus lowering the spectro- 
scopic binary detection efficiency. However we directly 
examine the stellar companion fraction as a function of 
spectral type or color B — V in Fig. 11. For both single- 
lined and double-lined spectroscopic binaries, if the bi- 
nary detection efficiency was systemically higher for K 



dwarfs then the anti-correlation could be a selection ef- 
fect. However we find that it is fairly independent of 
spectral type. Thus the anti-correlation does not appear 
to be a spectral- type dependent selection effect. 

We also examine the spectral type (~ mass) depen- 
dence of the correlation between planetary companions 
and host metallicity. The linear best-fit to the plan- 
etary companions of the FG sample has a gradient of 
a' = (0.22 ± 0.09)/5.7% = 3.9 ± 1.6 and the non-linear 
best-fit is a — 2.3 ± 0.6 as shown in Fig. 8 for the 
close d < 25 pc stars. These are significant at the 2 
and 3 a levels respectively. The linear best-fit to the 
planetary companions of the K sample has a gradient 
of a! ='(0.11 ± 0.10)/4.0% = 2.8 ± 2.5 and the non- 
linear best-fit is a. — 1.6 ± 0.9 as shown in Fig. 9 for the 
close d < 25 pc stars. These are both significant at be- 
tween the 1 and 2 a levels. The K sample contains fewer 
planetary and stellar companions compared to the FG 
sample. Both the linear and non-linear fits are consis- 
tent between the FG and K samples suggesting that the 
correlation between the presence of planetary compan- 
ions and host metallicity is independent of spectral type 
and consequently host mass. The fraction of planetary 
companions is also fairly independent of spectral type as 
shown in Fig. 11. 

Thus our results suggest that the correlation between 
the presence of planetary companions and host metal- 
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Fig. 10. — We compare the linear a' (triangles) and non-linear a (circles) parameterizations for the various samples listed in Table 3. The 
red points are the best-fits to planetary companions and the green points the best-fits to stellar companions. The fact that the red planet 
values for a are significantly larger than zero confirms and quantifies the metallicity/planet correlation. The fact that the green stellar 
values for a are predominantly less than zero, significantly so only for K dwarfs, is a surprising new result. The labels on the RHS refer to 
the samples for which the best-fit parameterizations are valid. We normalize the linear parametrization by dividing the best-fit gradient 
a by the average companion fraction Pavg (see text). This plot is a graphical version of Table 3 whose notes also apply to this plot. All 
of the best-fits are from [Fe/H] plots except for our FGK stars where a is the average of the best-fits to both the [Fe/H] and Z/Zq plots. 
The Pq values plotted in the vertical panel on the right refer to the corresponding best-fit normalization at solar metallicity (Eqs. 2-5). 



licity is significant at the ~ 4c level and that the anti- 
correlation between the presence of stellar companions 
and host metallicity is significant at the ~ 2a for the 
d < 25 pc FGK sample. Splitting both samples into FG 
and K spectral type stars suggests that the correlation 



between the presence of planetary companions and host 
metallicity is independent of spectral type but that the 
anti-correlation between the presence of stellar compan- 
ions and host metallicity is a strong function of spectral 
type with the anti-correlation disappearing for the bluer 
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Fig. 11. — Color (B — V) distribution for double-lined (squares) 
and single-lined (triangles) spectroscopic binaries (SB2s and SBIs 
respectively) and exoplancts (circles) in our close d < 25 pc sample. 
The linear best-fit gradient for SB2s is 0.00 ± 0.06, for SBIs it is 
-0.08±0.08 and for exoplanets it is -0.05±0.08. All three of these 
gradients are only significant at ^ la level. There is no significant 
correlation between SB1, SB2 or planetary fraction for either FG 
(B — V < 0.75) stars or for K (B - V > 0.75) stars. 



FG host stars (see Fig. 10). We find no spectral- type de- 
pendent binary detection efficiency bias that can explain 
this anti-correlation. 



4. IS THE ANTI-CORRELATION BETWEEN 
METALLICITY AND STELLAR BINARITY 
REAL? 

We further examine the relationship between stellar 
metallicity and binarity by comparing our sample with 
that of the Geneva and Copenhagen survey (GC) of the 
solar neighbourhood (Nordstrom et al. 2004) that has 
been selected as a magnitude-limited sample, a volume- 
limited portion (d < 40 pc) of which we analyse. This 
selection criteria infers that the sample is kinematically 
unbiased, i.e., the sample contains the same proportion 
of thin, thick and halo stars as is found in the solar neigh- 
bourhood. We also compare our sample with that of the 
"Carney-Latham" survey (CL) that has been kinemati- 
cally selected to have high proper motion stars (Carney 
& Latham 1987), i.e., it contains a larger proportion of 
halo stars compared to disk stars than is observed for the 
solar neighbourhood. 

Our sample is based on the Hipparcos sample that 
has a limiting magnitude for completeness of V = 
7.9 + 1.1 sin \b\ (Reid 2002) where b is Galactic latitude. 
Thus the Hipparcos sample is more complete for stars 
at higher Galactic latitudes where the proportion of halo 
stars to disk stars increases. Hence our more distant 
(25 < d < 50 pc) sample will have a small kinematic 
bias in that it will have an excess of halo stars, whereas 
our closer (d < 25 pc) sample will be less kinematically 
biased. 
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Fig. 12. — Histogram of stars in the complete volume limited GC 
sample (d < 40 pc). We exclude those stars with b — y< 0.3 so that 
the spectral type range becomes F7-K3 and thus similar to that of 
our sample. We only include those stars that have between 2 and 
10 radial velocity measurements with the CORAVEL spectrograph. 
We find an anti-correlation between binarity and host metallicity 
as shown by the linear and non-linear best-fits represented by the 
dashed and solid lines respectively. 



4.1. Comparison with a Kinematically Unbiased 
Sample 

The GC sample contains primarily F and G dwarfs 
with apparent visual magnitudes V 5; 9 and is complete 
in volume for d < 40 pc for F0-K3 spectral type stars. 
Unlike our sample analyzed in Section 3, it also includes 
early F spectral type stars. The GC sample color range 
is defined in terms oib—y not B — V like our samples. We 
remove these early F stars with b — y < 0.3 (B — V < 0.5, 
Cox 2000) from the sample so that the GC sample spec- 
tral type range is similar to ours. The GC sample then 
ranges from 0.3 < b - y < 0.6 (0.5 < B - V < 1.0) with 
those stars above b — y = 0.5 (B — V ~ 0.75) referred to 
as K stars. We also exclude suspected giants from the 
GC sample. 

For the GC sample, we only include those binaries ob- 
served by CORAVEL between 2 and 10 times so as to 
avoid a potential bias where low metallicity stars were 
observed more often, thus leading to a higher efficiency 
for finding binaries around these stars. This homoge- 
nizes the binary detection efficiency such that any real 
signal will not be removed by such a procedure. Un- 
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Fig. 13. — Same as Fig. 12 but only for the FG dwarfs in the 
GC sample of stars (d < 40 pc). We define FG dwarfs as those 
with b — y < 0.5 (B — V £ 0.75). We find only a marginal anti- 
correlation between binarity and host metallicity as shown by the 
linear best-fit with gradient a = —0.05 ± 0.06 and the non-linear 
best-fit with a = -0.11 ± 0.10. Using P avg = 26.1%, the scaled 
linear gradient a' = (-0.05 ± 0.06)/26.1% = -0.2 ± 0.2. 

like our sample for which we only include binaries with 
P < 5 years, the GC sample also includes much longer 
period visual binaries in addition to short period spectro- 
scopic binaries such that the total binary fraction of all 
types corresponds to ~ 25%. Comparing this with the 
period distribution for G dwarf stars of Duquennoy & 
Mayor (1991) this binary fraction corresponds to binary 
systems with periods less than ~ 10 5 days. 

For the volume limited d < 40 pc sample we again 
find an anti-correlation between binarity and stellar host 
metallicity as shown in Fig. 12. Both the linear and 
non-linear best-fits listed in Table 3 are significant at or 
above the 3a level. We also split the GC sample into FG 
and K spectral type stars in Figs. 13 and 14 respectively. 
The anti-correlation between the presence of stellar com- 
panions and host metallicity is significant at less than the 
la level for FG stars but significant at the ~ 4cr level for 
K stars. 

These results are qualitatively the same as those found 
for our sample but quantitatively weaker as shown in Fig. 
10 (confer rows of points labeled GC). This may be due 
to the higher fraction of late F and early G spectral type 
stars compared to our samples or the larger range (~ 
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Fig. 14. — Same as Fig. 12 but only for the K dwarfs in 
the GC sample of stars (d < 40 pc). We define K dwarfs as 
those with b — y > 0.5 (B — V S 0.75). We find a very strong 
anti-correlation between binarity and host metallicity as shown 
by the linear best-fit with gradient a = —0.46 ± 0.10. The non- 
linear best-fit is a = —0.52 ± 0.11. The scaled linear gradient 
a' = (-0.46 ± 0.10) /32.0% = -1.4±0.3. Comparing this plot with 
Fig. 13 suggests that the anti-correlation between binarity and 
host metallicity is stronger for redder stars. 

10 5 days) in binary periods contained in the GC sample 
compared to our sample where P < 5 years. Another 
way of interpreting this anti-correlation between binarity 
and metallicity may be in terms of the age and nature of 
different components of the Galaxy described by stellar 
kinematics, i.e., F stars are generally younger than K 
stars and thus are more likely to belong to the younger 
thin disk star population than the older thick disk star 
population. Hence we examine our results in terms of 
stellar kinematics. 

4.2. Comparison with a Kinematically Biased Sample 

We also compare our samples and that of the GC sur- 
vey with the Carney & Latham (1987) high proper mo- 
tion survey (CL). The CL survey contains all of the A, 
F and early G, many of the late G and some of the early 
K dwarfs from the Lowell Proper Motion Catalog (Gi- 
clas et al. 1971, 1978) and which were also contained in 
the NLTT Catalog (Luyten 1979, 1980). The number of 
stars in this distribution increases as the stellar colors be- 
come redder, peaking at about B — V = 0.65, following 
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which the numbers of stars begin to decrease (Carney 
et al. 1994). This group has also obtained data for a 
smaller number of stars from the sample of Ryan (1989) 
who sampled sub-dwarfs (metal-poor stars beneath the 
main-sequence) that have a high fraction of halo stars in 
the range 0.35 < B— V < 1.0. We refer to this combined 
sample as outlined in Carney et al. (2005) as the CL 
sample. This CL sample contains all binaries detected as 
spectroscopic binaries, visual binaries or common proper 
motion pairs. 

In Fig. 15 we plot the binary fraction of stars on pro- 
grade and retrograde Galactic orbits as shown in Fig. 
3 of Carney et al. (2005). All of the CL stars have 
[Fe/H] < 0.0. The CL distribution contains a small sub- 
set of metal-poor [Fe/H] < —0.2 stars from Ryan (1989) 
that has a one-third lower prograde binary fraction due 
to fewer observations. Thus stars with metallicities of 
between —0.2 and 0.0 have a higher binary fraction than 
the rest of the CL distribution. We make a small cor- 
rection for this bias in the binary fraction in the range 
—0.2 < [Fe/H] < 0.0 by lowering the 2 highest metallicity 
prograde points of the CL distribution by 2%. 

We note an anti-correlation between the binary frac- 
tion and metallicity for —1.3 < [Fe/H] < 0.0 range of 
prograde disk stars of the CL distribution as shown in 
Fig. 15. We find that the linear best-fit to this anti- 
correlation has a gradient of a = —0.07 ± 0.06 and the 
non-linear best-fit has a = —0.12 ± 0.09, which are both 
significant at slightly above the la level. For consistency 
we exclude the two lowest metallicity points from this 
best-fit so that we analyse the same region of metallicity 
as our samples and the GC sample and because these 
two low metallicity points will probably contain a signif- 
icant fraction of halo stars. The average binary fraction 
is Pavg = 26% for the disk-dominated part of the pro- 
grade CL distribution. Carney et al. (2005) found no 
correlation between binarity and host metallicity for the 
retrograde halo stars. 

We overplot our d < 25 pc binary fraction (from 
Fig. 4) along with the GC d < 40 pc binary frac- 
tion (from Fig. 12) onto the prograde CL sample in 
Fig. 15. All three of these samples have different bi- 
nary period ranges and levels of completeness. We scale 
our sample and the GC sample to the size of the Car- 
ney et al. (2005) sample by scaling the distributions 
to contain the same number of binary stars at solar 
metallicity. The most metal-poor point in our close bi- 
nary distribution is scaled above 100%, hence we set 
this point to 100%. The combined three sample distri- 
bution shows an anti-correlation between binarity and 
metallicity. The normalized linear best-fit to this is 
a' = (-0.10 ± 0.03)/25.7% = -0.39 ± 0.12 and the non- 
linear best-fit is a = —0.22 ± 0.05 which are both signifi- 
cant at or above the ~ 3cr level (see last row of Table 3). 
This combined result is our best estimate and indicates 
a strong anti-correlation between stellar companions and 
metallicity for [Fe/H] > -1.3. 

4.3. Discussion 

We examine our results in terms of Galactic popula- 
tions by determining the most likely population mem- 
bership (halo, thick or thin disk) for each star in the GC 
sample using the method outlined in the Appendix and 
then plotting them in the Galactic tangential velocity V 
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Fig. 15. — This plot is adapted from Fig. 3 of Carney et al. 
(2005). The black triangles are the points from the CL sample of 
proper motion stars with prograde Galactic tangential velocities. 
We overplot the binary fraction as a function of host metallicity 
for our close (d < 25 pc) F7-K3 sample (Fig. 4) with red circles 
and the green squares are from the volume limited GC sample 
(d < 40 pc) for F7-K3 stars (Fig. 12). The three samples contain 
different average binary fractions because the period range and the 
levels of completeness of the stellar companions varies between the 
samples as discussed in the text. We normalise the distributions 
by scaling our sample and the GC sample so that they contain the 
same fraction of binary stars as the sample of Carney et al. (2005) 
at [Fe/H] = 0. The linear and non-linear best-fits to the three 
samples combined are shown as dashed and solid lines respectively. 



- metallicity [Fe/H] plane as in Fig. 16. We use red 
points for the thin disk stars, green points for the thick 
disk stars and a blue point for the single halo star. The 
kinematically unbiased GC sample contains mostly thin 
disk stars. Excluding the one halo star in the GC sample, 
stars with [Fe/H] < — 0.9 belong to the thick disk and 
stars with [Fe/H] > — 0.1 belong to the thin disk. The 
region —0.9 ^ [Fc/H] > — 0.1 contains a combination of 
both thick and thin disk stars. 

We also plot these thick and thin disk stars as sepa- 
rate histograms in metallicity in Fig. 17. In the region 
[Fe/H] < — 0.9 that contains only thick disk stars, we 
find that the binary fraction is approximately twice as 
large as for the region [Fe/H] > — 0.1 that contains only 
thin disk stars. In both of these single population regions 
the binary fraction also appears to be approximately in- 
dependent of metallicity. While our purely probabilis- 
tic method of assigning the stars in the GC sample to 
Galactic populations is useful for determining the gen- 
eral regions of parameter space that the individual pop- 
ulations occupy, it is not precise enough to show exactly 
which stars belong to which population. This is espe- 
cially true for the regions of parameter space that have 
large overlaps such as that between the thick and thin 
disk stars in Fig. 16. Thus the thin and thick disk bi- 
nary fractions in the interval —0.9 ^ [Fe/H] < — 0.1 are 
probably mixtures. We suspect that the thin and thick 
disk binary fractions in this overlap region will remain at 
the same levels as found for the non-overlapping regions. 
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Fig. 16. — We plot tangential Galactic velocity V as a function 
of metallicity [Fe/H] for the kinematically unbiased GC sample 
(d < 40 pc). We use a probabilistic method to assign the stars in 
the GC sample to the three Galactic populations (halo, thick and 
thin disks) as discussed in the Appendix. Red points are thin disk 
stars, green points are thick disk stars and a blue point is the single 
halo star in the sample at V < —500 km/s. Crosses represent FG 
spectral- type stars and circles K stars. The ratio of thick/thin disk 
stars is ~ 3 times higher for K stars than for FG stars. 



The anti-correlation between binarity and metallicity in 
the —0.9 < [Fe/H] < — 0.1 range may be due to this over- 
lap between higher binarity thick disk stars and lower 
binarity thin disk stars. 

We now partition the Galactic tangential velocity V - 
metallicity [Fe/H] parameter space into four quadrants. 
We split the V parameter space into those stars on pro- 
grade Galactic orbits (P) and those on retrograde Galac- 
tic orbits (R). We split the [Fe/H] parameter space into 
those stars that are metal rich (r) with [Fe/H] > — 0.9 
and those that are metal poor (p) with [Fe/H] < — 0.9. 
We then label these quadrants by the direction of Galac- 
tic orbital motion followed by the range in metallicity or 
Pp, Pr, Rr and Rp as shown in Fig. 16. We now as- 
sume that the Pp quadrant contains a mixture of halo 
and thick disk stars and that the Pr quadrant contains 
a mixture of thin and thick disk stars and that the Rp 
and Rr quadrants only contain halo stars. 

The combined anti-correlation between binarity and 
metallicity in Fig. 15, that all three samples appear to 
have in common, is predominantly in the Pr quadrant of 
V — [Fc/H] parameter space that contains a mixture of 
thick and thin disk stars. As discussed above this anti- 
correlation may be due to the overlap of high binarity 
thick disk stars and lower binarity thin disk stars. 

While Latham et al. (2002) suggest that the halo and 
disk populations have the same binary fraction, Carney 
et al. (2005) find lower binarity in retrograde stars. As 
shown in Fig. 15 there is a clear difference of about a fac- 
tor of 2 in the region [Fe/H] > — 0.9 between the binary 
fractions of prograde disk stars compared to retrograde 
halo stars (Pr and Rr respectively). All the retrograde 
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Fig. 17. — Histogram of the stars in Fig. 16 suspected of belong- 
ing to the thick disk and the thin disk in the GC sample (d < 40 
pc). Notice the difference by a factor ~ 2 between the higher bi- 
nary fraction thick disk stars and the lower binary fraction thin disk 
stars. We note that the K star distribution contains a higher ratio 
of thick disk stars than the FG star distribution in the thick/thin 
disk overlap region. 



halo stars appear to have the same binary fraction (quad- 
rants Rr and Rp). The Pp quadrant contains prograde 
halo stars and has a ~ 2 times higher binary fraction than 
the quadrants containing retrograde halo stars. However 
the Pp quadrant also contains thick disk stars in addition 
to prograde halo stars. 

We propose that the Pp region, [Fe/H] < — 0.9, for 
stars on prograde Galactic orbits contains a mixture of 
low binarity halo stars and high binarity thick disk stars. 
In Fig. 15 at [Fe/H] ~ -0.9 our close sample (d < 25 
pc) and the GC sample (d < 40 pc) start to diverge 
from the data points of the CL survey. This observed 
divergence may be due to the CL survey being comprised 
of high proper motion stars and consequently a higher 
fraction of prograde halo stars compared to thick disk 
stars than the kinematically unbiased GC sample and 
our relatively kinematically unbiased sample where thick 
disk stars probably numerically dominate over halo stars. 

Using a kinemically unbiased sample, Chiba & Beers 
(2000) report that for the three regions —1.0 > [Fe/H] > 
-1.7, -1.7 > [Fe/H] > -2.2 and -2.2 > [Fe/H] that 
the fraction of stars that belong to the thick disk are 
29%, 8% and 5% respectively, with the rest belonging to 
the halo. We restrict these thick disk fraction estimates 
to stars only on prograde orbits by assuming that all 
of the thick disk stars are on prograde orbits and that 
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half of the halo stars are prograde and the other half 
are retrograde. Thus the fraction of prograde stars that 
are thick disk stars is 45%, 15% and 10% for the three 
metallicity regions respectively. 

Using these three prograde restricted thick disk/halo 
ratios reported in Chiba & Beers (2000) combined with 
observed binary fraction for the thick disk (55%) and 
halo (12%) stars in Fig. 15 we can test the proposal in 
the Pp quadrant that the two lowest metallicity prograde 
points from Carney et al. (2005) contain a mixture of low 
binarity halo stars and high binarity thick disk stars. We 
plot the three estimated mixed thick disk/halo binary 
fraction points as grey triangles in Fig. 15. We note 
that they are consistent with the two prograde Carney 
et al. (2005) points thus supporting our proposal that 
the Pp quadrant contains a mixture of low binarity halo 
stars and high binarity thick disk stars. These mixed 
thick disk/halo points also show a correlation between 
the presence of stellar companions and metallicity for 
stars in the Pp region. 

Our results suggest that thick disk stars have a higher 
binary fraction than thin disk stars which in turn have a 
higher binary fraction than halo stars. Thus for stars on 
prograde Galactic orbits we observe an anti-correlation 
between binarity and metallicity for the region of metal- 
licity [Fe/H] > —0.9 that contains an overlap between the 
lower-binarity, higher-metallicity thin disk stars and the 
highcr-binarity, lower-metallicity thick disk stars. We 
also find for stars on prograde Galactic orbits, a cor- 
relation between binarity and metallicity for the range 
[Fe/H] < — 0.9, that contains an overlap between the 
highcr-binarity, higher-metallicity thick disk stars and 
the lower-binarity, lower-metallicity halo stars. 

5. SUMMARY 

We examine the relationship between Sun-like (FGK 
dwarfs) host metallicity and the frequency of close com- 
panions (orbital period < 5 years) . We find a correlation 
at the ~ 4cr level between host metallicity and the pres- 
ence of planetary companion and an anti-correlation at 
the ~ 2a level between host metallicity and the presence 
of a stellar companion. We find that the non-linear best- 
fit is a = 2.2 ± 0.5 and a = —0.8 ± 0.4 for planetary and 
stellar companions respectively (see Table 3). 

Fischer & Valenti (2005) also quantify the planet 
metallicity correlation by fitting an exponential to a his- 
togram in [Fe/H]. They find a best-fit of a = 2.0. Our 
result of a = 2.2±0.3 is a slightly more positive correla- 
tion and is consistent with theirs. Our estimate is based 
on the average of the best-fits to the metallicity data 
binned both as a function of [Fc/H] and Z/Z . Larger 
bins tend to smooth out the steep turn up at high [Fe/H] 
and may be responsible for their estimate being slightly 
lower. 

We also analyze the sample of Nordstrom et al. (2004) 
and again find an anti-correlation between metallicity 
and close stellar companions for this larger period range. 
We also find that K dwarf host stars have a stronger anti- 
correlation between host metallicity and binarity than 
FG dwarf stars. 

We compare our analysis with that of Carney et al. 
(2005) and find an alternative explanation for their re- 
ported binary frequency dichotomy between stars on pro- 
grade Galactic orbits with [Fe/H] < compared to stars 



on retrograde Galactic orbits with [Fe/H] < 0. We pro- 
pose that the region, [Fe/H] < —0.9, for stars on prograde 
Galactic orbits contains a mixture of low binarity halo 
stars and high binarity thick disk stars. Thick disk stars 
appear to have a ~ 2 higher binary fraction compared to 
thin disk stars, which in turn have a <~ 2 higher binary 
fraction than halo stars. 

While the ratio of thick/thin disk stars is ~ 3 times 
higher for K stars than for FG stars we only observe a 
marginal difference in their distributions as a function of 
metallicity. In the region —0.9^ [Fe/H] < — 0.1 that we 
suspect contains a mixture of thick and thin disk stars, 
the K star distribution contains a higher ratio of thick 
disk stars compared to the FG star distribution at a given 
metallicity. This difference is marginal but can partially 
explain the kinematic and spectral-type (~ mass) results. 

Thus for stars on prograde Galactic orbits as we move 
from low metallicity to high metallicity we move through 
low binarity halo stars to high binarity thick disk stars to 
medium binarity thin disk stars. Since halo, thick disk 
and thin disk stars are not discrete populations in metal- 
licity and contain considerable overlap, as we go from 
low metallicity to high metallicity for prograde stars, we 
firstly observe a correlation between binarity and metal- 
licity for the overlapping halo and thick disk stars and 
then an anti-correlation between binarity and metallicity 
for the overlapping thick and thin disk stars. 

We would like to thank Johan Holmberg for his help 
on analysing the Geneva sample and Chris Flynn, John 
Norris, Virginia Trimble, Richard Larson, Pavel Kroupa 
and David Latham for helpful discussions. 
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TABLE Al 

Properties of the Three Stellar Populations 



Component 




(V) 


ov 


a w 


<[Fe/H]> 


CT [Fe/H] 


Fraction / 


Thin Disc 


43 


-15 


28 


17 


-0.1 


0.2 


0.925 


Thick Disc 


67 


-53 


51 


42 


-0.8 


0.3 


0.070 


Halo 


131 


-226 


106 


85 


-1.8 


0.5 


0.005 



Note. — Data from Robin et al. (2003). 



APPENDIX 

PROBABILITY OF GALACTIC POPULATION MEMBERSHIP 

We use a similar method to that of Reddy et al. (2006) in assigning a probability to each star of being a member of 
the thin disk, thick disk or halo populations. We assume the GC sample is a mixture of the three populations. These 
populations are assumed to be represented by a Gaussian distribution for each of the 3 Galactic velocities U, V, W 
and for the mctallicity [Fe/H]. The age dependence of the quantities for the thin disk are ignored. The equations 
establishing the probability that a star belongs to the thin disk (Pthin), the thick disk (Pthick) or the halo (Phalo) arc 



Pthu 



fl p ) Pthick 



where 



/2-p ! Phalo 



f P 3 
73-p 



(Al) 



Pi = Ci exp 



1 



(A2) 



U 2 (V-(V)) 2 W 2 ([Fe/H] — ([Fe/H])) 2 



2 < 



ZO [Fc/H], 



< J u,&v i o'w i o'iF e /n] i 



-(i = l,2,3) 



Using the data in Table Al taken from Robin et al. (2003) we compute the probabilities for stars in the GC sample. 
For each star, we assign it to the population (thin disk, thick disk or halo) that has the highest probability. We plot 
the probable halo, thick and thin disk stars of the GC sample in Fig. 16. 



REFERENCES 



Amnions, S.M., Robinson, S.E., Strader, J., Laughlin, C, Fischer, 

D. & Wolf, A., 2006, ApJ, 638:1004-1017 
Batten, A.H., 1973, Pergamon, Oxford 

Bond, J.C., Tinney, C.G., Butler, P., Jones, H.R.A., Marcy, G.W., 

Penny, A.J. & Carter, B.D., 2006, MNRAS, 370:163-173 
Carney, B.W. & Latham, D.W., 1987, AJ 93:116-156 
Carney, B.W., Latham, D.W., Laird, J.B. & Aguilar, L.A., 1994, 

AJ 107:2240-2289 
Carney, B.W., Angular, L.A., Latham, D.W. & Laird, J.B., 2005, 

that Omega Centauri is Related to the Effect', ApJ, 129:1886- 

1905 

Cayrel de Strobel, C, Soubiran, O, Ralite, N., 2001, A&A, 
373:159-163 

Chiba, M. & Beers, T.C., 2000, ApJ, 119:2843-2865 
Cox, A.N., 2000, AIP Press, 4th Edition 

Dall, T.H., Bruntt, H. & Strassmeier, K.G., 2005, A&A, 444:573- 
583 

Duquennoy, A. & Mayor, M., 1991, A&A, 248:485-524 
Fischer, D.A. & Valenti, J., 2005, ApJ, 622:1102-1117 
Giclas, H.L., Burnham, R. & Thomas, N.G., 1971, Lowell 

Observatory, Flagstaff 
Giclas, H.L., Burnham, R., Jr. & Thomas, N.G., 1978, Lowell 

Observatory Bulletin, 8:89 
Gonzalez, G., 1997, MNRAS, 285:403-412 

Gonzalez, G., Laws, C, Tyagi, S. & Reddy, B.E., 2001, AJ, 121:432- 
452 

Grether, D. & Lineweaver, C.H., 2006, ApJ, 640:1051-1062 

Gunn, J.E. & Griffin, R.F., 1979, AJ, 84:752-773 

Hauck, B. & Mermilliod, M., 1998, A&A, 129:431-433 

Jones, H.R.A., Butler, P., Marcy, G.W., Tinney, C.G., Penny, A.J., 

McCarthy, C. & Carter, B.D., 2002, MNRAS, 337:1170-1178 
Jones, H., Butler, P., Tinney, C.H., Marcy, G., Carter, B., Penny, 

A., McCarthy, C.H. & Bailey, J. , 2006, MNRAS, 369:249-256 



Kotoneva, E., Flynn, C. & Jimenez, R., 2002, MNRAS, 335, 1147- 
1157 

Laughlin, G. & Adams, F.C., 1997, ApJ, 491:L51-L54 

Latham, D.W., Mazeh, T., Carney, B.W., McCrosky, R.E., 
Stefanik, R.P. & Davis, R.J., 1988, AJ, 96:567-587 

Latham, D.W., Stefanik, R.P., Torres, G., Davis, R.J., Mazeh, T., 
Carney, B.W., Laird, J.B. & Morse, J.A., 2002, ApJ, 124:1144- 
1161 

Latham, D.W., 2004, ASP Conference Series Vol 318, eds. Hilditch, 

R.W., Hensberge, H. & Pavlovski, K. 
Laws, C, Gonzalez, G., Walker, K.M., Tyagi, S., Dodsworth, J., 

Snider, K. & Suntzeff, N.B., 2003, AJ, 125:2664-2677 
Lineweaver, C.H. & Grether, D., 2003, ApJ, 598:1350-1360 
Luyten, W.J., 1979, University of Minnesota, Minneapolis 
Luyten, W.J., 1980, University of Minnesota, Minneapolis 
Martcll, S. & Smith, G.H., 2004, PASP, 116:920-925 
Nordstrom, B., Mayor, M., Anderson, J., Holmberg, J., Pont, F., 

j0rgensen, B.R., Olsen, E.H., Udry, S. & Mowlavi, N., 2004, 

A&A, 418:989-1019 
Reddy, B.E., Lambert, D.L. & Allendc Pricto, C, 2006, MNRAS, 

367:1329-1366 
Reid, I.N., 2002, PASP, 144:306-329 

Robin, A.C., Reyle, C, Derriere, S. & Picaud, S., 2003, A&A, 

409:523-540 
Ryan, S.G., 1989, AJ 98:1693-1767 

Santos, N.C., Israelian, Mayor, M., Bcnto, J.P., Almeida, P.C., 

Sousa, S.G. & Ecuvillon, A., 2005, A&A, 437:1127-1133 
Santos, N.C., Israelian, G. & Mayor, M., 2004, A&A, 415:1153-1166 
Valenti, J.A. & Fischer, D.A., 2005, ApJ, 159:141-166 



